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The treatment of long-range interacting systems (including Newtonian self- 
gravitating systems) remains a challenging issue in statistical mechanics. Due 
to the lack of extensivity, they present non-standard effects like negative specific 
heats, which shows the inequivalence of statistical ensembles (namely, microcanon- 
ical and canonical) even in the limit of infinite number of particles (TV — * oo). In 
this paper we review a series of results obtained for one and two dimensional simple 
TV-body dynamical models with infinite-range attractive interactions and without 
short distance singularities. The free energy of both models can be exactly obtained 
in the canonical ensemble, while information on the microcanonical ensemble and 
on the dynamical evolution can be derived from direct numerical simulations with 
simple O(TV) codes, which make use of mean-field variables. Both models show a 
phase transition from a low energy clustered phase to a high energy gaseous state, 
in analogy with the models introduced in the early 70's by Thirring and Hertel. 
The phase transition is second order for the ID model, first order for the 2D model. 
Negative specific heat appears in both models near the phase transition point, but 
while for the 2D model it is an equilibrium phenomenon, in the ID case it is typi- 
cal of transient metastable states, whose lifetime grows with TV. For both models, 
in the presence of a negative specific heat, a cluster of collapsed particles coex- 
ists with a halo of higher energy particles which perform long correlated flights, 
which lead to anomalous diffusion: the mean square displacement grows faster 
than linear with time. The dynamical origin of this " superdiffusion" is however 
different in the two models, being related to particle trapping and untrapping in 
the cluster in ID, while in 2D the channelling of particles in an egg-crate effective 
potential is responsible of the effect. Both models are Lyapunov unstable and the 
maximal Lyapunov exponent A has a peak just in the region preceeding the phase 
transition. Moreover, in the low energy limit A increases proportionally to the 
square root of the internal energy, while in the high energy region it vanishes as 
TV" 1 / 3 . Since the ID model is explicitely constructed considering the first modes 
of a Fourier expansion of a classical one dimensional gravity potential, the large 
scale properties of this model in the low-energy clustered phase resemble those of 
ID gravity (mass-sheet models). The relation of both models with gravity remains 
to be analysed in detail by adding successive Fourier waves. However, we believe 
that the main dynamical effects we have found: superdiffusion and "strong" chaos 
(due to the evaporation and sticking of particles to the cluster) should be present 
also in "true" gravitational systems, and we see no obstruction that they appear 
also in 3D. 

1 Introduction 

It is well known that extensivity is an essential ingredient for building up 
thermodynamics from statistical mechanics questions (see e.g. Landau's text- 



when internal energy, entropy and other thermodynamic quantities are no 
more extensive, i.e. when a part of a system has not the same thermody- 
namics properties of the whole. This problem naturally arises in Newtonian 
gravity, where the inherent long-range of the forces is the ultimate cause of 
non-extensivity, but is also present in plasma physics although the coexis- 
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about what happens 



tence of attractive and repulsive long-range Coulomb interactions determines 
the screening effect which mitigates this problem. Therefore, one knows that 
statistical physics does not have the same strong impact for gravitational 
systems as it has for systems with short-range interactions and hard-cores, 
although one would really like to use statistical concepts in this field^|. 

On the other hand mean-field models are frequently used in statistical 
mechanics as a first simple tool to understand the collective behavior of sys- 
tems with short-range interactions and it is believed that they describe their 
"universal" features (such as critical exponents at phase transitions) above a 
certain "critical" number of neighbors. In this context, it is curious that none 
really posed the question of the validity of the statistical approach, considering 
that mean-field models (expecially those written in terms of an infinite-range 
Hamiltonian) explicitly violate extensivity. A reason for that lies in the fact 
that the thermodynamic limit is performed resorting to saddle-point tech- 
niques, which put the Hamiltonian in the explicitly decoupled form, hiding 
the difficulties inherent in the long-range interaction. However, all concepts re- 
lated to statistical ensemble equivalence (like microcanonical/canonical equiv- 
alence) are open for such systems, and we will review in this paper models 
of this type for which it is indeed violated, producing effects like negative 
specific heats within the jpicrocanonical ensemble, which are well known in 
the gravitational context u. 

The central idea contained in the papers which we are reviewing here is 
that it is possible to establish a close correspondence between the behaviors 
observed in the Hamiltonian dynamics of iV-body self-gravitating systems 
and those of much simpler models of the mean-field infinite-range class. We 
have dealt only with systems in one and two spatial dimensions, but this 
approach can in principle be extended also to three dimensions. One of these 
models has been indeed rigorously derived from the self-gravitating model 
in one dimension, by restricting the line to a finite segment and truncating 
the Fourier expansion of the potential to its first wave component The 
two-dimensional model is less tightly related to two dimensional gravity, but 
arises from a similar truncation of the Fourier expansion in a box B. The 
truncation is not ineffective and introduces features which are not present in 
the original self-gravitating system. In one dimension, the truncated model 
displays a second order phase transition from a phase where the mass points 
are grouped in a cluster to one where they are homogeneously dispersed; 
this phase transition is not present in the full self-gravitating system, which 

a Recently, Tsallis and co-workers have developed a new interesting approach based on an 
alternative definition of entropy, which could cope with non-extensivityEI. The relation of 
their approach with our results remains to be carefully investigated 
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is always in the clustered phase (being the potential confining). However, 
we claim that many features of the clustered phase of our mean-field model 
are similar to those observed in the full model. Moreover, our model can 
be a toy system for studying phase transitions in the simple setting of one 
dimension (with all the computational advantages), since it is expected that 
phase transitions are present in three dimensions In two dimensional models 
some theoretical result U seem to suggest that a phase transition occurs at 
finite temperature from a clustered to a homogeneous phase. 

The advantage of studying one and two-dimensional mean-field models 
is that their free energy in the canonical ensemble can be exactly derived by 
performing the mean- field limit (the infinite N limit at fixed volume), which 
is a reasonable surrogate of the thermodynamic limit (the infinite -jY n limit 
at fixed density). Unfortunately, at variance with Thirring models &□, the 
microcanonical solution (for instance the entropy-energy relation) is not yet 
available for such models (although some attempts already exist along this 
line 113). However, plenty of informations can be gotten from relatively fast 
numerical simulations, thanks to the algorithmic advantage of working with 
O(N) codes (instead of usual 0(N 2 )), due to the introduction of mean-field 
variables. Moreover, important informations can be obtained from solving 
the one-dimensional non collisional Boltzmann-Poisson (BP) equation Ej, al- 
though it corresponds to inverting the N — > oo limit with the infinite time 
limit, which can be a dangerous exchange in the presence of transient states. 
!_ . In particular, from the analysis of the BP equation for the ID model 
113, it turns out that all the initial conditions with a homogeneous particle 
distribution and with a symmetric velocity distribution, are stable in time. 
Moreover, extensive dynamical simulations of the Id and 2D model allowed us 
to study in detail aspects like non-ergodicity and long-time memory, that axe 
expected phenomena in systems with non-screened long-range interactions t3. 

The structure of this paper is the following. After introducing and dis- 
cussing the main models in Section 2, we derive statistical equilibrium prop- 
erties in Section 3. In this Section we comment also on the presence of 
metastable states near the phase transition point in the one-dimensional 
model; these states have a lifetime which grows with N. In the two- 
dimensional model, the microcanonical simulations reveal a region of negative 
specific heat. Dynamical properties are discussed in Section 4. The main 
effect is anomalous diffusion, superdiffusion in our case, which means that the 
mean square displacement of a single point mass grows in time faster than lin- 
ear. In the one-dimensional model this phenomenon appears as a consequence 
of the presence of metastable states, i.e. it is a transient-to-equilibrium effect 
and is generated by the correlated trapping and untrapping of the particles in 
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the cluster. In the two-dimensional model it is instead due to the generation 
of an effective mean-field potential of the egg-crate form, where particles in 
a certain energy range can perform long free flights; the phenomenon is de- 
stroyed by noise on the long-time span, but noise dyes out as N — > oo and thus 
superdiffusion extends to the time asymptotic regime. In Section 5 it is shown 
that the dynamics of both models is chaotic, by computing the maximal Lya- 
punov exponents A. Since the systems are integrable for vanishing internal 
energy U and in the high energy limit, A should exhibit a peak at intermedi- 
ate energy. Indeed, A is maximal in the phase transition region. In the high 
energy regime A vanishes with an universal scaling law A = X(N) ~ AT -1 / 3 , 
while at low energy it is positive also in the mean-field limit and increases as 
A = \(U) ~ [7 1 / 2 . In Section 6 we finally draw some conclusions and discuss 
some perspectives of this novel approach. 

2 The Models 

In this Section we introduce the ID and 2D models. We discuss their Hamil- 
tonians and equations of motion by introducing the mean-field variables. We 
also give a brief description of the phenomenology and of the tools that have 
been introduced for their study. 

2.1 One dimensional model 

This first model describes a collection of N identical particles that move 
under their mutual gravitational-like attraction. Each particle is a mass point 
and boundary conditions are periodic. Particles are therefore constrained to 
move on the unit circle, z_e. in a one dimensional geometry. The total energy 
is constant and given by □: 



where 9% 6] — tt, n] is the coordinate of the i-particle, pi its conjugate momen- 
tum and m the mass of the particles. A time-discrete version of this model 
has been previously introduced by Kaneko and Konishi lia. The system is 
isolated, total momentum is conserved and can be set to zero without loss of 
generality. K and V are the kinetic and potential energy respectively. As the 
latter has no singularity when i = j, the particles do not collide but smoothly 
cross through one other. This shows up clearly in the equations of motion 
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(1) 



N 

= £ = (2) 

since they correspond to a system of fully coupled pendula. 

The potential energy is rescaled by 1/N in order to get a finite specific 
energy TL = H/N in the thermodynamic limit N — > oo, hence the potential is 
thcrmodynamically stable (H > EqN). The characteristic period t c of parti- 
cle motion is then an intensive quantity, given by t c = 27r/y / ro. In the large 
N limit, particles can be treated in a mean-field approximation and the cUstci- 
bution functions can be computed using Poisson-Boltzmann techniques t3t5. 
These techniques have been used in a large variety of dynamical systems (plas- 
mas, galaxies, fluids) where particle-particle collisions can be neglected when 
compared to the particle-mean-field coupling. For model (|l|), the mean-field 
variable writes: 



M = (M x , My) = j^^Sj ' with S J = (cos0,-,sin0j) . (3) 
With these definitions, Eq. (0) reduces to that of a single perturbed pendulum 



§i = -Msm(6i-<j>) (4) 

where we have set the mass to = 1 for simplicity and = arg(M) is the phase 
of the mean-field. Both M and <f> are scalar quantities that depend on time 
through 9' s and their time evolution describes the collective dynamics of the 
system. 

The single-particle dynamics is ruled by the Hamiltonian 

2 

tH = ^ + l- M x cos(e t ) - My sm(9i ) = K l + V l (5) 

which is non authonomous, since it depends on the mean-field, which reflects 
the collective dynamics of the model. Hence, each particle moves in a mean- 
field potential Vj determined by the instantaneous positions of all the other 
particles of the system. 

The thermodynamics of model (^) can be derived theoretically in the 
canonical ensemble (see Sect. 3). It shows a second-order phase transition at 
a critical value of the internal energy U c — H c /N = 0.75 between a collapsed 
phase, where a single cluster is present, and a homogeneous phase, where 
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particles are uniformly distributed on the circle. The development of the 
cluster in the collapsed phase-corresponds to the well known Jeans instability 
for self-gravitating systems til 

Numerical simulations of self-gravitating systems in their collapsed phase 
and initially far from equilibrium, usually reveal two evolution mechanisms: 
i) collective phenomena that involve almost all the particles and that rapidly 
drive the system into a quasi-stationary configuration; ii) individual effects 
that rule the successive slow relaxation to equilibrium. On short time scales, 
spiraL structures are often developed to eliminate the excess of kinetic en- 
ergy □; they are due to the formation phase-space domains with an excess or 
depletion of particles. These structures are robust enough to resist individual 
particle effects and are known to prevent the system from reaching equilib- 
rium They are similar to the Dupree structures observed in plasmas where 
they inhibiL_Landau damping and constitute an important element of Vlasov 
turbulence liH. 

Model (0) contains only one spatial scale, the length 2tt; this reflects in the 
presence of a single modulation of the particle density (the cluster). However, 
this model belongs to the family of TV-body Hamiltonian systems that write 



where k n = 
harmonics I 



1=1 i,j=l n=Q n \ J 

+- 1 and where the parameter s counts the number of Fourier 
The interest in this family of Hamiltonians lies on the fact 



that in the limit s — > oo they write 



N 2 2 N N 2 2 

where G is the gravitational constant and the customary dimensionless units 
1 = 2nG have been chosen Ej. Up to the 1 /N scaling constant of the potential, 
this Hamiltonian is the one of a collection of N identical planar mass sheets of 
constant mass surface density /x evolving in one dimension perpendicular to 
their surface. This model was originally introduced to describe the dynamics 
of stars evolving in the neighborhood of a flattened galaxy El For the Poisson- 
Boltmann limit to exist, one must require the the total mass fiN — m remains 
finite when one lets N — > oo. Then 

TT N -2 

i—1 i T j>i 
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where pi — Pi/N are the scaled momenta. We would like to notice that 
the Hamiltonian can be transformed in that corresponding to the usual 
gravitational ID model also by simply multipling the time units by a factor 
uV, this obviously corresponds also to rescale the energy units by a factor N 



The dynamics of the mass sheet model involves an interplay of all the 
spatial scales, which is the source of its metastability. The advantage of us- 
ing models (|^) is that they allow a selection of a limited number of modes, 
through the tuning of the parameter s; the qualitative changes in the dynam- 
ical and in the thermodynamical behavior can be followed by adding higher 
Fourier components. Such a multiscale approach is based on the mean-fields 
Mfcn = s j,K with Sj t k n = (cos k n 6j, sin k n 8j). This approach to self- 

gravitating systems is similar to the one proposed in Ref. E3 where a weak 
formulation based on a smoothed-out gravitational potential is introduced. 
We expect that the large scale properties won't be strongly modified by the 
introduction of smaller spatial scales. However, the balance between collective 
(large scale) and individual (small scale) effects might depend crucially on the 
number of Fourier components. The investigation of the properties of model 
(|l|) is thus the very first step in the investigation of Hamiltonian (ft). Indeed, 
one major difference is that while model ([!]) has a second-order phase transi- 
tions, model (@) has no phase transition. The absence of phase transitions is 
a trivial consequence of the existence of the scale transformation 9i —> a 2 9i, 
Pi — > api, t — > t/a, which transforms the energy as H.^ —> c^Tioo and leaves 
the equations of motion unchanged. Thus, if a > 1 states corresponding to 
a larger energy coincide with those at a shorter time. Therefore, increasing 
the energy cannot produce any phase transition, but it is like observing the 
dynamics on a shorter time scale. Hamiltonian (j^) is always in the clustered 
phase, since the potential is always confining. To induce a phase transition 



one must explicitly introduce a scale in the model (this 
for the spherically concentric mass shell model in Ref. 
is how the phase transition disappears in the s — > oo limit, and a related oncis 
whether phase transitions are present in 2D and 3D self-gravitating systems la. 



,s been recently done 
). An open question 



2.2 Two dimensional model 

We have also studied a two dimensional model, that can be viewed as a gener- 
alization of the one dimensional model introduced in the previous Subsection. 
We consider a system of N pointlike identical particles with unitary mass 
moving in a 2-D box with periodic boundary coditions. The dynamics is 
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ruled by the Hamiltonian 

N v 2 +» 2 

H = l^ 9 



(=1 

1 N 
-Y 



2N 



3 - cos(a;i - Xj) - cos(yi - yj) - cos(x l - Xj) cos(j/i - j/j) 



= K + V , 



(9) 



where (x^yt) e] - 7r,7r]x] - 7r,7r], (xi,p i}X ) and (yi,p i}V ) are pairs of conju- 
gate variables, if (resp. V) is the kinetic (resp. potential) energy. The 
reference energy is chosen in order to ensure the vanishing of the total en- 
ergy of the system when all the particles are placed in the same position 
(V = 0) at rest (K = 0). Without the presence of the third term in the po- 
tential the particles would move independently along the x and y-directions 
according to the previously introduced 1-D potential (|l|). The interparticle 
potential appearing in ([)]) mimicks the Fourier expansion in a periodic square 
box of side L = 2ir of a 2-D self-gravitating Newtonian potential (log |r|, be- 
ing r the interparticle Euclidean distance), restricted to its first three terms 
(n x ,n y ) — (0, 1), (1, 0), (1, 1), where k x = 2im x /L and k y = 2im y /L are the 
wave numbers along the two directions. Interactions of the log \r\ type arises 
also for point vortices in 2D Eulerian turbulence c3, and our approach could 
be fruitfully extended to this case. 

The equations of motion for the coordinates 



1 N 

iV ^ 



3 = 1 



sm(xi — xj) + sin(xi — Xj) cos(yi — yj) 



(10) 



and those for yi are obtained exchanging x «-> y, due to the symmetry of (|9|). 

Previous investigations of model ^ a have revealed that at low energy 
U = H/N the particles, due to attractive coupling, collapse into a unique 
cluster, with [/-dependent spatial extension. This clustered phase survives 
up to a critical energy U c ~ 2.0 (numerically determined from an implicit 
equation derived in the canonical ensemble), where a first-order phase tran- 
sition occurs to a homogeneous phase. This transition will be discussed in 
more detail in the next Section. The collective behaviour is evidenced by the 
following mean-field variables 

Mi = (< cos(a;) > N , < shi(x) > N ) = M x (cos(^i), sin(0i)) (11) 



M 2 = (< cos(y) >at, < sin(y) > N ) = M 2 (cos(^ 2 ),sin(0 2 )) (12) 



paperoJ^: submitted to World Scientific on February 1, 2008 9 



Pi = (< cos(x + y) > N , < sin(x + y) > N ) = P x (cos(-0i), sin(V'i)) (13) 



P 2 = (< cos{x - y) > N , < sin(x - y) > N ) = P 2 (cos(-0 2 ), sin(V> 2 )) (14) 

where < .. >n denotes the average over all the particles in the system. The 
moduli Mi,2 and Pi j2 are maximal and equal to 1 when all the particles have 
the same position and their value decrease when the spatial distribution of the 
particles extends. For a homogeneous distribution, according to the central 
limit theorem, Mi, 2 w Px,2 ~ 0(l/y/~N). These quantities are the order 
parameters characterizing the degree of clustering of the system. 

By reexpressing the equation of motion ( |l0| ) for both coordinates x and 
y in terms of the mean-field variables, one straighforwardly shows that the 
time evolution of the i— th particle is ruled in a way similar to the ID model 
by the following single particle non autonomous Hamiltonian : 



Since V is invariant under the transformations x <-> —x, y <-* — y and x <-> y, 
it turns out that in the mean-field limit (N — > oo with constant U = H/N), 
Mi = M 2 = M and P 1 = P 2 = P. 

3 Equilibrium properties 

In this Section we will mainly discuss the canonical equilibrium exact solutions 
of the previoulsy introduced ID ([!]) and 2D (|^) self-gravitating models and we 
will compare these solutions with the result of the microcanonical simulations. 
The canonical free energy is easily obtained in both cases by resorting to an in- 
verse Gaussian transformation (sometimes called the Hubbard-Stratonovich 
trick) and evaluating an integral with saddle-point techniques. For the ID 
model (Q), we have shown that the phase transition between the collapsed 
and the homogeneous phase is a second order phase transition. For the 2D 
model (||) the transition is instead of first order. Microcanonical simulations 
reveal the presence of a negative specific heat region near the phase transition 
for both the ID and the 2D model, but for the former case this is only a tran- 
sient effect (however, the transient grows with the number N of point masses). 
Ensemble inequivalence and negative specific heats are indeed rather familiar 
to astrophysicts and appear in the phenomenon known as the gravothermal 
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= x ' 1 v ' % + 3 - Mi cos(xi - 4>i) - M 2 cos(yi - <t> 2 ) 

-ipi cos(a; l + yi - tpi) - ^P 2 cos(xi -yi-tpz) 
= Ki + Vi . 




(15) 



catastrophe These observations confirm the conjecture of Lynden-Bell Eat 
and demontrate on a simple model that the equivalence among statistical 
ensembles is violated due to the long-range nature of the forces and to the 
absence of repulsive hard-cores. 

3.1 One dimensional model 

The one-dimensional model of Eq. ([!]) can be exactly solved in the canonical 
ensemble. Its free energy 

F = - lim —In lYlde.dpiexpi-H/T) (16) 

i 

is given by the expression 

- F/T = ~ \n(2nT) - |1 + niax^^ + \n(2irI Q {x))} (17) 
where the max is obtained by solving the consistency equation 

£4 = ^, (18) 

where Iq and I\ are respectively the zero and first order modified Bessel 
functions. In the following we will adopt the value m = 1, without loss of 



generality. Eq. (18) has the unique solution x — 0, corresponding to the 
vanishing mean-field region M = 0, for T > 0.5; it has instead two symmetric 
nonvanishing solutions, corresponding to a nonvanishing value of M for T < 
0.5. This latter is the clustered phase and the ratio M/T determines the 
average size of the cluster (which becomes infinitely narrow, i.e. it collapses 
to a point, only in the zero temperature limit). The temperature dependence 
of the mean field M is exactly derivable and it is shown in the inset of Fig. [l] 
by the continuous line. The temperature-energy relation can be also derived 
using the standard formula U(T) = d(F/T)/d(\/T) and is reported in the 
same figure. This model has therefore a second order phase transition at the 
critical temperature T c = 0.5, corresponding to the critical energy U c = 3/4. 
The dynamical time averages of the same quantities can be determined by 
computer simulations. Being the model Hamiltonian, the energy is fixed by 
the initial state and temperature is dynamically determined by computing the 
average kinetic energy per particle, and waiting for long time relaxation. 

In all energy ranges, except close to the critical point, every chosen initial 
state relaxes quite fast to the canonical curves of M vs. U and T vs. U ; 
finite N corrections are visible in Fig. [l]. On the contrary, in a region just 
below the critical point (evidenciated by the vertical line in Fig. |l|) initial 
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Figure 1. Temperature T versus the internal energy U for the ID model (|l|). The full curve 
is obtained theoretically in the canonical ensemble (TV — > oo). The points are the result of 
microcanonical simulations with an increasing number of mass points, in particular the full 
circles refer to an initial "water bag distribution". In the inset, the mean- field variable M 
vs. U is reported with the same notation. The horizontal lines indicates U c - 



states of the "water-bag" class (uniformly distributed momenta in a fixed 
interval symmetric around zero and all positions at 9i = 0) show a relaxation 
to metastable non-equilibrium states, whose T vs. U curve has two main 
features: i) for U > 0.67 the points lye on the continuation for U < 3/4 of the 
line T = 2U — 1; the corresponding states are gas- like, with zero magnetization 
(uniform in space) and have non-Gaussian tails in momentum; ii) for 0.6 < 
U < 0.67, T decreases when U increases, this is the negative specific heat 
region. All these latter points, obtained for a system of N = 20, 000 point 
masses, relax very slowly to the canonical curve (as shown by the behavior of 
their relative Boltzmann entropy in Fig. |^a) in a time which increases linearly 
with the number N of particles of the system Fig. |^b. Ideally, if one would 
have prepared a system with an infinite number of particles, this would not 
have relaxed to the canonical equilibrium. The metastability of our model 
is not as strong as the one encountered for disordered systems (where the 
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relaxation time diverges exponentially in TV) , it is of the same type as the one 
found for the mass-sheet model and for plasma models as reviewed in Ref. t3. 
Its origin lies in the complex structure of the chaotic web through which the 
orbit channels in phase-space (see further comments Section 4.1). 



0.20 




0e+00 2e+03 4e+03 6e+03 8e+03 1e+04 1e+04 

N 



Figure 2. a) Absolute difference with respect to the equilibrium Boltzmann entropy for the 
momentum distribution vs. time t. Relaxation slows down for increasing N. b) A plot of 
the relaxation time t p vs. the number of particles N. The full line is a best linear fit to the 
data. 

The equilibrium behavior discussed above and depicted in Fig. ^ can be 
derived also using the non-collisional Boltzmann-Poisson (BP) equation in two 
different ways: i) imposing the maximal entropy principle, as done in Ref. EH; 
ii) looking for solutions with a Gaussian shape in momentum and a factorized 
stationary distribution function f(9,p,t) — /i(#)/2(p) ta- To our knowledge, 
the fact that the canonical equilibrium solution coincides with the one given 
by the BP equation is a novel feature, which deserves further investigations 
for other models. The use of the BP equation gives access also to the one- 
particle distribution function f{6,p,t), on which numerical checks have been 
performed in RefEj. A feature which remains unexplored in the context of the 
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PB equation is metastability, but recently one of the authors t3 has .proven, 
using the expression of the evolution operator in terms of Lie brackets 1^3, that 
initial states symmetric in momentum fi(p) = fi(—p) and uniform in position 
do not evolve in time to all orders in perturbation theory. This is a first sign 
of metastability in the PB context, where since one is treating directly the 
N — » oo limit, metastability manifests itself in a multi-stability. Finally, let 
us mention that recently in Ref. EH it has been shown that the transition 
observed for the model ([!]) is associated, in the mean- field limit, to a change 
in the topology of the corresponding configuration space. 

3.2 Two dimensional model 

Model (^) can be solved in the canonical ensemble as in the ID and all its 
equilibrium features can-be exactly derived in the mean- field limit, using again 
saddle-point technique u. We do not reporte here tha-explicit expression for 
the free energy, that is already reported elsewhere The single particle 
potential V% for the 2D case (|l5|) is reported in Fig-.||, it has a shape similar to 
the so-called egg-crate potential studied in Ref. C3. The potential is periodic 
along the two spatial directions and in each elementary cell has 4 maxima 
(Ym — 3 + 2M — P), 4 saddle points (V s = 3 + P) and a minimum (V m = 
3 — 2M — P). The depth of the potential well is (V s — V m ) and the center of 
the cluster in the collapsed phase coincide with the position of the minimum 
of the potential. In the limit U — ► 0, M and P — ► 1 and therefore Vm — > 4, 
V s — > 3 and V m — > 0. In this limiting situation all the particles are trapped in 
the self-consistent potential well. Increasing U, the kinetic energy produces 
an evaporation of the particles from the cluster. This implies that also the 
values of M, P and of the well depth decrease (see Fig. ||). For U > U c , there 
are no more clustered particles and in the limit N — > oo the single particle 
potential becomes flat Vm — V s — V m — 3 and time independent. 

However, due to finite N effects, the instantaneous mean-field variables 
fluctuate with a typical time 01 within the statistical band range ~ l/ViV. 
This implies that in the collapsed phase the particles move in a fluctuating 
potential. 

In Fig.^the temperature T is reported as a function of U. Above U c ~ 2.0, 
T increases linearly with U indicating that the system behaves like a free 
particle gas. In the collapsed phase, the tendencsz_jof the system to collapse 
is balanced by the increase of the kinetic energy cB. This competition leads 
initially (for < U < 1.8) to a steady increase of T, followed (for 1.8 < U < 
U c ) by a rapid decay of T. This yields a negative specific heat as illustrated 
in the inset of FigJ|. These results are in full agreement with theoretical 
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Figure 3. Instantaneous single particle potential V(x,y) together with its contour-plot for 
U = 1.00 and (a;, y) £ [-5tt/2, 5tt/2] 2 . 



predictions based cp-the analysis of a simple classical cell model a and with 
numerical findings ESEJ for short ranged attractive potentials without hard- 
cores. The phenomenon of negative specific heat can be explained within 
a microcanonical approach with an heuristic argument u: approaching the 
transition, a small increase of U leads to a significative reduction of the number 
of collapsed particles (as confirmed from the drop exhibited by M and P for 
U > 1.8 in Figji]); as a consequence the value of V grows and, due to energy 
conservation, the system becomes cooler. 

nn Our data confirm also another important prediction of Hertel and Thirring 
an: the non-equivalence of canonical and microcanonical ensemble nearby 
the transition region. In the inset of Fig. we report the microcanonical 
findings, obtained by direct simulations, and the theoretical canonical results. 
These two sets of data coincide everywhere, except in the energy interval 
1.6 < U < 2.0. The discrepancy is due to the impossibility of the canonical 
ensemble to exhibit a negative specific heat, a prohibition that does not hold 
for the microcanonical ensemble. Our theoretical estimation of the Helmholtz 
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Figure 4. Meanfield quantities M and P as a function of U. The solid line refer to canonical 
results analytically estimated, while the symbols to microcanonical results obtained via 
numerical simulations. The exponent a defined in eq. (Ed), are also shown (asteriskes). 
The data have been obtained with N = 4, 000 (apart few points with TV = 10, 000) and 
averaged over a total integration time of order 10 6 . 



free energy F — F(T) reveals that usually F has an unique minimum. For 
T < 0.5 the minimum F c corresponds to non zero values of M and P (i.e. 
to the collapsed phase), while for T > 0.55 the minimum Fjj is associated to 
M = P = (i.e. to the homogeneous phase). In the region 0.5 < T < 0.55, 
both minima F c and Fjj coexist as local minima of the free energy. However, 
for T < T c = 0.54 the collapsed phase is observed because F c < Fji , while 
for T > T c the homogeneous phase prevails since Fjj < F c . At T = T c 
the two minima are equivalent and a jump in energy from U(T~~Lst 1.6 to 
£/(T+) « 2.0 is observed. This transition is therefore first order E3, with a 
finite latent heat U(T+) - U{T~). 

At variance with the ID model, there is no sign of metastability in the 
2D (apart eventually at U c ), at least for what concerns the thermodynamical 
properties. All the data reported in Figs. ^ and || correspond to the time 
asymptotic state, also in the negative specific heat region. This is the main 
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u 

Figure 5. Temperature T as a function of the internal energy U for the 2D model. The 
solid line refers to analytical results obtained within the canonical ensemble, while the full 
circles to microcanonical simulation results. The inset is an enlargement of the transition 
region: the full curves (resp. dashed) refer to the absolute (resp. relative) minimum of the 
free energy F(T). 



not yet understood difference between the ID and the 2D model. 
4 Dynamical Properties 

In this Section we discuss the main dynamical properties of the ID and 2D 
models. The most important phenomenon is correlated to trapping and un- 
trapping of particles within the cluster in the low energy collapsed phase, this 
effect produces an anomalous diffusion (superdiffusion) of particles, which has 
a finite time span, being destroyed by relaxation to equilibrium in the ID 
model and by noisy fluctuactions (due to the finite number of particles) of the 
single-particle potential in the 2D model. However, for both models in the 
mean-field limit superdiffusion should be present at any time. 

4-1 One dimensional model 

The single particle Hamiltonian (5) of the ID model is that of a perturbed 
pendulum. The perturbation has two different sources in the time dependence 
of M, which generates the stochastic layer around the separatrix, and in that 



papero4: submitted to World Scientific on February 1, 2008 17 



of cj), which is the main origin of the trapping/untrapping phenomenon in the 
nonlinear resonance of the pendulum. In Fig. ^ we show a snapshot taken 
at some large time of the particles which were initially put in a "water-bag" 
initial condition: near the center of the resonance particles librate forming a 
spiral whose branches touch the separatrix region. The unperturbed separa- 
trix is drawn making reference to the asymptotic value of M and its center 
in momentum is slightly displaced from zero, because the cluster is slowly 
drifting □. Out-of-resonance particles are capable of performing long flights if 
energy is large enough. A typical trajectory of a particle is shown in Fig. ^ 
and reveals the presence of tracts of almost free motion with a typical positive 
or negative speed, interrupted by tracts where the particle is trapped in t 
cluster, performing with it a much slower motion. It has been shown in Ref. 
that the motion is " superdiffusive" , i.e. the mean square displacement < 2 > 
(where the average is taken over the particles and over distinct time origins) 
has the anomalous scaling 

< e 2 >~ t a (19) 

with a — 1.4±0.1 in a given energy range just below the phase transition point, 
0.6 < U < 0.7, and for a finite time. When the system relaxes to thermal 
equilibrium, diffusion becomes normal (a = 1), but the time taken to do this 
change increases with N, exactly as the relaxation time to equilibrium. We 
can then state that in the thermodynamic limit the motion is superdiffusive 
forever, when starting from a non-equilibrium (e.g. "water-bag") initial state. 
It is suggestive, and deserves further investigation, the fact that superdiffusive 
motion is found in the region of negative specific heat, where a core of collapsed 
trapped particles coexists with a "halo" of almost free particles. 



4-2 Two dimensional model 

We concentrate in this Subsection on the description of the transport proper- 
ties of the 2D model (|J). In the collapsed phase each particle of our system 
movesJn the single particle potential Vi = Vi(t) (see Fig. ^) of the egg-crate 
form c3. For a single particle moving in such fixed potential the self dynam- 
ics is known to be anomalously diffusive when the particle is channeling, i.e. 
when its energy lies between Vm and V s E3. This anomalous behaviour is 
due to the competition of laminar and localized phases. Indeed, channeling 
particles intermittently show an almost ballistic motion along the channels of 
the potential interrupted by localized sequences, where the particles bounces 
back and forth on the maxima of the potential. 

Usually anomalous diffusion has been studied in systems with a few de- 
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Figure 6. Snapshot at some large time of the single-particle phase-space for an initial "water 
bag" distribution at U = 0.4 and with N = 10, 000 for the ID model. The full line indicates 
the separatrix of the unperturbed pendulum drawn at the corresponding equilibrium value 
of M. 



grees of freedom b±3t£3. Qritvijew attemps have been made to consider TV-body 
dynamics with N ^> 1 oE2a. 

In this Subsection we give some indications relative to the basic dynamical 
mechanisms governing single particle transport. To this aim, we consider the 
time dependence of the mean square displacement (MSQD), that usually reads 
as 

<r 2 (t)>oct a (20) 

where the average < . > is performed over different time origins and over 
all the particles of the system. The transport is said to be anomalous when 



a =/= I: namely, it is s 
and ballistic for a = 2 



lse if < a < 1, superdiffusivc if 1 < a < 2 
1. The usual Einstein diffusion law corresponds 
to a = 1 and in 2D can be written as < r 2 (i) >= 4Dt, where D is the 
self-diffusion coefficient. 

If one considers a single particle with an initial energy between Vm and 
V s and follows its trajectory for some time it displays features quite similar to 
the so-called Levy walks Ej (see Fig.||). This kind of trajectories are usually 
identified when anomalous diffusion occurs. 



papero4: submitted to World Scientific on February 1, 2008 19 




Figure 7. Free flights and trapped motions of a particle in the fD model for U = 0.6, 
N = 500 and a "water bag" initial condition. 



As we already mentioned, finite iV-effects play a relevant role in the dy- 
namics. Due to self-consistency, they are responsible for the fluctuations in 
time of the mean-field quantities Mi. 2 and Pi. 2- The potential experienced 
by each particle thus fluctuates in time and particles having an energy close 
to V s have the possibility to be trapped in the potential as well as to escape 
from it. This implies that the localization phenomena illustrated in Fig are 
not only due to bounces of the particle on the maxima of the potential, but 
also to trapping in the potential well itself due to separatrix crossing. 

We have argued from simple considerations that diffusion should be 
anomalous in the collapsed phase. This is indeed the case, as confirmed from 
direct evaluations of the MSQD in a quite extended energy interval. An ex- 
ample of this is reported in Fig.^| for U — 1.1 and N = 4,000. Diffusion is 
anomalous for times smaller than a crossover time r beyond which the Ein- 
stein's diffusion law is recovered < r 2 (t) >= ADt. A similar behaviour for the 
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Figure 8. ^Typical time evolution of a particle initially in a channel of the 2D single particle 
potential ha. In this representation, the 2 — D torus on which the dynamics takes place is 
unfolded and represented as an infinite plane constituted of an infinite number of juxtaposed 
elementary periodic cell of size lxl. In the inset an enlargement of the trajectory in the 
indicated box is reported. 



MSQD was previously .observed for a system of N symplectic (globally and 
locally) coupled maps but with a subdiffusive (i.e. with a < 1) short 

time dynamics. 

The energy dependence of the a-values is illustrated in Fig. [|. It shows 
up clearly that the thermodynamical phase transition from collapsed to ho- 
mogeneous phase is associated to a dynamical transition from superdiffusion 
(with 1.3 < a < 1.9 for 0.4 < U < 2.0) to ballistic motion (with a ~ 2 for 
U > U c ). In the collapsed regime, we observe an increase of a from 1.3 ± 0.1 
to 1.9 ± 0.1, that is due to the modification of the shape of the single particle 
potential. 

A possible explanation of this behaviour is the following. At very low en- 
ergy U < 0.3, all the particles are essentially trapped in the potential well and 
none is channeling: no diffusion is observed in this case. For increasing energy 
a fraction of particles (due to the decrease of Vm — V m ) escapes from the well 
and some of them get enough energy to move along the channels: anomalous 
diffusion is then evidenced. The increase in the value of the exponent a is 
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Figure 9. Time dependence of the mean square displacement < r 2 (t) > in a log-log scale 
for U = 1.1 and N = 4, 000 for the 2D model. The solid line indicate the simulation results. 
The dot-dashed segments are the estimated slopes of < r 2 (t) > for t < r and t > r. 



due to the fact that also the channels width Vm — V s grows with U. However, 
for U approaching U c the number of untrapped particles increases noticeably, 
but now the channel width vanishes abrubtly. This implies that a significant 
fraction of particles will move freely (with energy > Vm)- These mechanisms 
lead naturally to ballistic motion for U > U c , where the potential Vi is now 
almost constant apart from fluctuations of order 0(l/y/~N). 

The fact that in the asymptotic limit (t — > oo) normaLdiffusion is recov- 
ered constitutes a typical signature of a noisy dynamics E3 ! E3H. But in our 
system no external noise is added to the system, therefore the transition from 
anomalous to asymptotic ordinary diffusion must be attributed to a "deter- 
ministic source of noise" , that is intrinsic in our system and due to finite size 
effects. 

A deterministic anomalously diffusing dynamical system submitted to 
weak environmental white noise shows a transition from anomalous diffusion 
to standard diffusion on sufficiently long time scales E3. The crossover time 
r for this transition can be explicitely computed l£I and turns out to increase 
as an inverse power law of the noise amplitude, when short time behaviour is 
superdiffusive. 
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When analyzing the N dependence of r for two energies, U = 1.48 and 
U = 2.0, we find that r oc N. The interpretation of the N dependence 
of the crossover time is straighforward if we consider finite N effects as a 
source of noise in our model, with a typical amplitude of order 1/yiV. This 
last assumption is justified by the fact that the microscopic dynamics of the 
particles generates stochastic fluctuations 0(1/VN) in the values of M and 
P. These fluctuations become weaker for increasing N and thus naturally 
yield an increasing value of r. Mope details about the relevant dynamical 
mechanisms can be found elsewhere u. 

We would like to stress that in the mean-field limit the diffusion will 
be anomalous at any time, while for finite N normal diffusion will be always 
recovered asymptotically. The fact that the dynamics of the system is strongly 
influenced by the order in which the two limits N — > oo and t — > oo are taken 
should be related to the long-range nature of the forces present in our model. 

5 Chaotic properties 

In order to complete the description of the two models, we investigate in this 
Section a fundamental indicator that characterizes the dynamics of Hamil- 
tonian systems: the maximal Lyapunov exponent A. Our analysis relies on 
numerical estimations of A, performed considering the evolution in thc-.tangent 
space and applying the standard technique first introduced in Ref. E3. 

Our models are integrable in the limit of low and high energy, therefore 
A — > both for U — » and U — > oo. In between these two extrema we expect 
that a finite Lyapunov exponent will be observed. 

5.1 One dimensional model 

The ID model is Lyapunov unstable (A > 0) over all the energy range, see 
Fig. [19 (the maximal Lyapunov exponent of model (^) was first computed in 
Ref. Below U ~ 0.2 there is a weak size dependence of the Lyapunov 
exponent and a substantial independence on the chosen initial condition. A 
scaling law A ~ J7 1 / 2 is clearly visible (see Fig. fLl|a) and still lacks a convincing 
theoretical explanation (the theoretical expression derived in Ref. 113 gives 
A ~ U,la sharp disagreement with numerical data); an argument, reported 
in Ref. Ii3 heuristically justifies the correct scaling. For U > 0.2, A rapidly 
increases and reveals a strong dependence on N. In fact, it can be shown that 
in all the vanishing M phase (U > U c ), A vanishes as iV -1 / 3 (see Fig. pT|b)). 
This scaling is easily obtained resorting to a random matrix approximation 
of the dynamics 113, which considers the successive Jacobians as random and 
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statistically independent, and has been also derived theoretically in Ref. 
using differential geometric techniques. Moreover, also in true one dimensional 
gravitational models A vanishes as an inverse power law of N as shown in EjO 
(in those cases the exponent ranges from 1/5 to 1/4). For true 3D gravitational 
models Gurzadyan and Savvidy, by relating the average Riemaniann curvature 
to the corresponding Lyapunov spectrum, were able to suggest the following 
scaling law c2l : 

GMN 2 / 3 , N 

AG « ^ (21) 

As already reported in Sec. 2.1, our models can be connected to true gravi- 
tational models by simply changing the time scale by a factor 1 / VN, since in 
our model the mass M is unitary, the gravitational constant G is dimension- 
less, and the volume V — (27r) d (where d is the dimension of the system) do 
not depend on N. Therefore it is sufficient in Eq.([2l|) to rescale the average 
velocity of the particles < v > and the Lyapunov exponent by a factor y/N. 
After such rescaling one obtains from Eq.(pl|) A cx TV -1 / 3 . This suggests the 
validity of our results also in the 3D gravitational case. The maximal Lya- 
punov exponent relaxes to its asympotic value with a N~ 1 ^ 3 -\aw also for a 
gas of hard-spheres lB: this universality would deserve further investigation. 

Although the strong N dependence hides the effect, an evident peak of 
A is present just below the phase transition point, in correspondence of the 
negative specific heat region. Chaos is "stronger" in this region due to the con- 
tribution of particles which evaporate from the cluster, entering the stochastic 
layer near the separatrix. The curious coexistence of negative specific heat, 
superdiffusion and maximal chaos deserves further studies. A confirmation of 
such a behavior of the maximal Lyapunpy exponent for systems with long- 
range interactions can be found in Ref. Ell, where, slightly modifying model 
(|]), it has been shown that when the interaction becomes short-ranged A no 
more vanishes in the N — + oo limit. These authors further argue that the 
exponent —1/3 of the infinite-range model modifies continuously by reduc- 
ing the range of the interaction, reaching zero when the interaction becomes 
short-ranged. 

5.2 Two dimensional model 



Our data for the 2D model (g) are reported in Fig. [12| for three different types 
of initial conditions : 

• (A) Maxwellian velocity distribution and clustered particles; 
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Figure ID. Maximal Lyapunov exponent A as a function of the internal energy for the ID 
model (|l|) for increasing N. The system has been started from equilibrated initial data. 



• (B) Maxwellian velocity distribution, but with a thermal velocity coin- 
ciding with its canonical prediction, and particles spatially organized in 
a single cluster in such a way that also the M- and P-values reproduce 
the canonical prediction; 

• (C) water-bag velocity distribution and clustered particles. 

The initial condition (C) is the one commonly used through the present 
paper, in particular for the study of transport in the system. As it can be 
seen from Fig. ^| , A grows for increasing U up to a maximum value and 
then decreases. Such maximum (at least for N = 200) is located at an energy 
U ~ 1.3 — 1.4 < U c . For U > U c , we observe a power law decrease of A with 
N with an exponent ~ 0.31 in good agreement with the results for the ID 
model. 
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Figure 11. a) Maximal Lyapunov exponent A vs. the internal energy U for the ID model 
various systems sizes N, showing TV-independence and the f/ 1 / 2 common scaling, b) A 
vs. N for various values of U > U c = 0.75. For increasing value sof U the scaling law 
jy— 1/3 j s ver jfi ec [ better and better. The stars refer to a simulation performed with random 
Jacobians, showing the effectiveness of the random matrix approximation. 



In our model the vanishing of A in the mean-field limit, for U > U c , is 
connected to the fact that the single particle potential becomes flat and the 
system integrable. 

In the low energy limit (for U < 0.01) a power law increase of the type 
U 1 / 2 is clearly observable for all the three types of initial conditions similarly 
to ID. This indicates that for fully coupled Hamiltonian systems this property 
holds in general, independently of the space dimensionality. In particular, for 
U — > the particles are all clustered, therefore we expect that the scaling 
A oc C/ 1 / 2 should be related to a "collective" chaotic mechanism. 

Let us now discuss the mechanisms behind the observed behaviours of 
A. Once the energy U is fixed for all the three considered initial conditions, 
after a reasonable transient, we obtain exactly the same values for the aver- 
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Figure 12. Maximal Lyapunov exponent as a function of U in a log-linear plot, for three 
different initializations of the 2D model. The dashed line indicates the critical energy U c . 
The reported data corresponds to N = 200 and to integration times ranging from t ~ 1 X 10 6 
to t ~ 9 X 10 6 . 



age temperature T and magnetizations M and P and a common Maxwellian 
distribution for the velocities. However, at low energies (0.001 < U < 0.8) the 
measured average A depends heavily on the initial conditions. This clearly in- 
dicates the coexistence of several equivalent states, that can be considered as 
equilibrated within the examined time interval. Evidently the measurement 
of A is more sensible to the presence of such equivalent states than that of 
thcrmodynamical variables. It should be noticed that this kind of behaviour is 
unexpected in iV-body Hamiltonian systems, because it is commonly believed 
that for sufficiently large values of N Arnold diffusion takes place and each 
orbit is allowed to visit the full phase-space. But our data instead indicate 
that some "barrier" in the phase space still survive even for N = 200. The 
origin of this lack of ergodicity is related to the long range nature of the forces 
that induces a persistent memory, of the initial conditions, as previously no- 
ticed by Prigogine and Severne 113 for gravitational plasmas. Recently, some 
numerical evidence of non ergodicity has been reported also for ID gravity 
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1 1 But ina recent paper El this evidence is questioned in favour of a slow 
relaxation mechanism affecting high energy particles. 
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Figure 13. Time evolution of the maximal Lyapunov exponent A, of the magnetization M, 
and of the kinetic energy K are shown for a initial condition of (B)-type for the 2D model 
at U = 0.87 and with N = 200. 

It is clear from Fig.|l2| that, in the interval U € [0.001,0.8], for initial 
conditions of type (B) A remains always smaller than the corresponding expo- 
nents obtained with initial conditions (A) and (C). The maximal differences 
are observed in the energy range 0.2 < U < 0.8, where particles begin to es- 
cape from the cluster (this happens for initial conditions (A) and (C)). Above 
U ~ 0.9 the same Lyapunov is obtained for all types of initial conditions. A 
typical feature of initial conditions (B), for U < 0.8 is that all the particles are 
trapped in the potential well. Instead when one or more particle escape from 
the cluster (U > 0.9) also with this initial condition the common value of A is 
obtained. Indeed, two chaotic mechanisms are present in the system: one felt 
by the particles moving in the minimum of the potential and one by particles 
visiting a region near to the separatrix. This second mechanism is well known 
and it is related to the chaotic layer around the separatrix. The former one 
is due to the erratic motion of the minimum of the potential well. In order to 
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clearly identify such mechanisms we have followed the trajectory of a system 
initially prepared in condition (B) for an energy U — 0.87 and N — 200. On 
a short time scale all the particles are trapped and we measure an average 
value A ~ 0.13. When at a later time one particle escapes from the cluster A 
almost doubles its initial value (sec Fig. |l3|). In Fig. [l^ the magnetization M 
and the kinetic energy K are also shown. When the particle escapes from the 
cluster M shows a clear decrease as well as K . This last effect is due to the 
fact that potential energy V is minimal when all the particles are trapped, 
therefore if a particle escapes V increases and, due to energy conservation, 
K decreases. This is the phenomenon underlying the negative specific heat 
effect. From the above arguments we can identify a "strong" chaos felt from 
the particles approaching the separatrix and a "weak" chaos associated to the 
orbits trapped in the minimum of the potential. The presence of these two 
chaotic mechanisms, together with the non-ergodicity of the system, explains 
the strong dependence of the values of A on the initial conditions. 

6 Conclusions 

This paper reviews, discussing differences and similarities, a class of models 
of point masses without hard-cores and interacting with an infinite-range at- 
tractive potential in one and two dimensions. Both models may be seen as toy 
systems for studying long-range attractive interactions in an extremely simpli- 
fied setting. We consider them a dynamical extension of th&Jattice gas model 
introduced and studied by Hertel and Thirring 0. Thirring Em pointed out that 
a transition from a clustered to a homogeneous phase can be associated with 
a negative specific heat (within the microcanonical ensemble) either (I) when 
the potentials are thermodynamically unstable @ or (II) when the thermody- 
namical potentials are extensive quantities, but the interaction is long-ranged 
(at least in the limit TV — > oo). In a true gravitational potential both the 
above conditions are fulfilled. Studies,-af-iV-body models with short-range 
thermodynamically unstable potentials E3'E2l confirmed that condition (I) was 
sufficient to observe transition from a clustered to a gaseous phase, negative 
specific heat, inequivalence of canonical and microcanonical ensemble. Our 
studies concerns a continuous time model with thermodynamically stable po- 
tential (ff/v > EqN), but with long-range attractive interaction. Therefore 
they confirm that the characteristic associated to the transition observed by 
Hertel and Thirring for a simple lattice gas model are present also in more 
realistic models and that condition (II) is indeed a sufficient one. Let us try 
to summarize condition (I) and (II) in an unique prescription : in order to 
observe a transition belonging to the Hertel and Thirring universality class it 
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is necessary that the thermodynamical potentials (e.g. the entropy) avs-not 
additive even in the limit N — > oo and that the interaction is attractive e3. In 
other words if a system composed by an infinite number of particles is split- 
ted in two sub-systems typically any thermodynamical potential of the whole 
system coincide with the sum of the entropies of the two sub-systems, if this 
is not the case violation of the standard thermodynamics should be expected. 

Moreover, we show here that extremely simple infinite-range mean-field 
models display several interesting features, which are present also in gravity 
in one and two dimensions (we believe that our approach can be extended to 
three dimensions without conceptual difficulties). A drawback of our models 
is that the microcanonical solution has not yet been obtained exactly and 
we must rely on numerical simulations (apart from some results which can 
be obtained by the collisionless Boltzmann-Poisson equation) . This drawback 
turns however into an advantage if one is interested in dynamical effects, about 
which thermodynamics gives no hint. The main dynamical effect discussed 
in this review is " superdiffusion" (the particles mean square displacement 
grows faster than linear with time). Superdiffusion is present in the region 
of negative specific heat and is due to the particles which evaporate from the 
cluster and perform long flights before being trapped again (the configuration 
space is a torus). Evaporation is also the source of "strong" chaos: in the 
negative specific heat region the maximal Lyapunov exponent A is much larger 
than everywhere else. Although gravity has no phase transition in ID (the 
gravitational system is always clustered) and the question is open whether 
a phase transition is present in 2D 0, both the effects discussed above could 
be present in the negative specific heat region of self-gravitating systems, 
and we do not see any obstruction to their existence also in 3D. Finally, let 
us comment on a few scaling laws that we have found to be universal for 
all the models of the class we have introduced. In the gaseous phase the 
maximal Lyapunov exponent A vanishes as A -1 / 3 in analogy with what is 
found forpSnme specific gravitational systems and for other dynamical models 
of fluids \BB. In the low energy clustered phase A ~ U 1 / 2 . This should be 
true also in 3D and could be checked in "true" self-gravitating systems in 
their collapsed low-energy phase. 
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